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ABSTRACT 

We present new limits on the ejection of metal-rich old-population hyperve- 
locity stars from the Galactic center (GC) as probed by the SEGUE-2 survey. 
Our limits are a factor of 3-10 more stringent than previously reported, depend- 
ing on stellar type. Compared to the known population of B-star ejectees, there 
can be no more than 30 times more metal-rich old-population F/G stars ejected 
from the GC. Because B stars comprise a tiny fraction of a normal stellar pop- 
ulation, this places significant limits on a combination of the GC mass function 
and the ejection mechanism for hypervelocity stars. In the presence of a normal 
GC mass function, our results require an ejection mechanism that is about 5.5 
times more efficient at ejecting B-stars compared to low-mass F/G stars. 

1. Introduction 

Hypervelocity stars (HVS) have emerged as a promising way to probe the dynamics 
and physical conditions at the Galactic center (GC). Thus far, the discoveries of HVS have 
predominantly been B-stars, with the more recent addition of a small number of A-stars 
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f lBrown et all BOOa Edelmann et all boOa iHirsch et all bOQa Brown et al.l l2007aU bl. 120091 ). 
These objects have been detected at large distances from the GC (~ 70kpc); their flight 
times and spectral types are consistent with ages of less than 100-200 Myr. Empirically, the 
known HVS constitute a relatively young population. While the youth of this population 
may reflect a top-heavy mass function at the GC or an ejection mechanism that greatly 
prefers stars of this mass range, it is also possible that the increased difficulty of locating 
old-population0 HVS among a predominantly old halo population has distorted our picture of 
the conditions at the GC. In the absence of a complete kinematic census of the Galactic halo, 
it is substantially more difficult to find old-population HVS compared to young-population 
HVS. This is because the Galactic halo is composed of primarily old-population stars whose 
red colors are very similar to those of low-mass HVS, which make photometric pre-selection 
for spectroscopy non-trivial (Kollmeier & Gould 2007). By contrast, there are very few 
blue stars in the halo, just a handful that are blue due to their youth (such as "runaway B 
stars" ) and not many more old stars (like blue-horizontal branch stars) that live in the halo 
"legitimately" . Early type stars are more luminous and can be seen to larger distances, and 
thereby probe larger volumes, relative to the more typical G-stars that predominate in a 
normal mass function. These factors combine to make the photometric background of non- 
HVS to HVS significantly reduced for early-type stars. However, regardless of the practical 
difficulties, it is essential to find the old-population HVS (or determine that they do not exist) 
if one is to understand the physics at the GC (Kollmeier & Gould 2007). The power of the 
wide-field spectroscopic capability of the 2nd Sloan Extension for Galactic Understanding 
and Exploration (SEGUE-2) can be directly brought to bear on this problem. 

In a previous paper, we used the entire available Sloan Digital Sky Survey (SDSS) stellar 
database to place lim its on the ejection of metal-rich old-population hypervelocity stars 
( iKollmeier et al.ll2009l . hereafter K09). Using turnoff stars collected in SDSS throughout the 
lifetime of the survey, we analyzed the high-velocity, \Vq\ > 300 km s -1 , distribution of stars. 
The underlying sample contained nearly 300,000 stars — a motley collection of calibration 
stars, failed quasar targets, and fiber-fill. The sheer number of targets allowed us to obtain 
upper limits to the total ejection rate of T F < 60 Myr -1 and T G < 300 Myr -1 for unbound 
F and G sta r s, resp ectively. Comparing to estimates of the B-star ejection rate from the 



Brown et al.l ( I2007bl ) survey, we obtained a relative ejection rate of old to young metal-rich 



HVS consistent with a normal stellar mass function. To probe more deeply into this putative 
population, we have implemented a well-defined target selection algorithm within SEGUE-2 
to more cleanly search for metal-rich old-population HVS that have been ejected from the 



1 We use the term "old" as a short-hand for "long-lived". Similarly, "young" is used to denote "short- 
lived" stars. 
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GC. 

In Section 2, we describe SDSS, SEGUE, SEGUE-2, and our target selection. In Sec- 
tion 3 we present the results of this selection and derive our sensitivity to metal-rich old- 
population HVS from which our new limits for the ejection of this population are derived. 
In Section 4 we discuss these limits. 

2. Sample 

The sample we analyze here comprises a subset of objects targeted within the SEGUE-2 
survey. 



2.1. SDSS, SEGUE and SEGUE-2 



SDSS-I was an imaging and spectros copic survey that began routine operations in April 



2000 , and continued through June 2005 ( Fukugita et al 



20001 ; [Hogg et al.ll200ll; 



20041 : iGunn et al. 



2.5m telescope (IGunn et al. 



2006 



Lupton et al.ll2001 



1996 



Smith et al.ll2002t IStoughton et al 



Gunn et al. 



1998 



York et al. 



2002 



Ivezic et al 



Tucker et al.l 120061 ) . The SDSS and its extensions used a dedicated 



20061 ) located at the Apache Point Observatory in New Mexico. 



The Sloan Extension for Galactic Understanding and Exploration (SEGUE) is one of the 
three key projects in the recently completed first extension of the Sloan Digital Sky Survey, 
known collectively as SDSS-II. The SEGUE program, which ran from July 2005 to July 2008, 
obtained ugriz imaging of approximately 3500 deg 2 of sky outside of the SDSS-I footprint, 
with special attention being given to scans of lower Galactic latitudes (|6| < 35°) in order 
to better probe the disk/halo interface of the Milky Way. SEGUE obtained approximately 
240,000 medium-resolution spectra of Gal actic stars, selecte d to explore the nature of stel- 
lar populations from 0.5 kpc to 100 kpc flYanny et al.l 120091 ). SDSS-III, which is presently 
underway, has already completed the sub-survey SEGUE-2, an extension primarily intended 
to obtain stellar spectra for distant stars that are likely to be members of the outer Galactic 
halo. 



The first seven public data releases from SDSS ( lAbazajian et al.ll2003l 12004 120051 12009 



Adelman-McCarthy et al.ll2006l 120071 . 120081 ) have produced over 350,000 stellar spectra (and 
their derived atmospheric parameters, where possible). More than 120,000 stellar spectra 
obtained during the course of SEGUE-2 will be distributed as part of the next public data 
release, DR8. 



The SEGUE Stellar Parameter Pipeline processes the wavelength- and flux-calibrated 
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spect ra generated by the standard SDSS spectroscopic reduction pipeline (jStoughton et al. 



20021 ). obtains equivalent widths and/or line indices for more than 80 atomic and molecular 



absorption lines, and estimates radial velocities, TW, log g , and [Fe/H] thro ugh the applica- 
tion of a number of approaches ( lAUende-Prieto et al.ll2008l ; iLee et al.ll2008al Jbl). The spectral 
resolution of the survey is R~2000 over wavelengths 385-920 nm, yielding typical velocity 
errors of 4-10 kms -1 . 



2.2. Target Selection 



While the target selection method advocated by iKollmeier &: Gould I (120071 ) is relatively 
complete, it would require many fibers on each SEGUE-2 plate. To achieve greater econ- 
omy of fibers, we adopted the following rather stringent selection for HVS candidates in 
SEGUE-2. The new selection criteria make use of proper motion information as well as 
photometric metallicity information to pre-select stars that have enhanced likelihood of be- 
ing GC ejectees. This experiment, de signed to lo cate the "metal-rich contaminant" to the 
Galactic halo originally envisioned by iHillsl (119881 ) . selects stars to have a high probability 
of being high-metallicity, and to have proper motions consistent with high space velocities 
from the GC. We select stars satisfying the following criteria: 



0) \b\ > 30° 

1) 17 < g < 20 

II) fi > 8 masyr" 1 

III) v tot > 400 kms" 1 

IV) either 

- IVa) < 6masyr _1 or 

- IVb) v± > 400 kms" 1 

V) either 

- Va) (0.35 < (g - r) < 0.40) n (0.375 <{u- g) - 2.5(g - r) < 0.525) or 

- Vb) (0.40 <{g- r) < 0.60) n (0.225 <{u- g) - 2.5(g - r) < 0.425) 



where n is the total pro per motion as determined by SDSS and USNO-B (IGould &: Kollmeier 
20041 ; iMunn et al.ll2004l ). fi± is the component of proper motion perpendicular to the direction 
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from the GC to the star, v to t is the total (3D) velocity (assuming that the star is on a radial 
orbit), and v± = fid is the transverse velocity assuming that the measured proper motion is 
correct (velocities quoted here are Galactocentric). There are several other criteria that are 
applied to all SEGUE-2 targets involving proper motions that are designed to remove poorly 
measured objects due to errors in USNO-B. These are described in greater detail below. 
However, the crucial point at this stage is that these criteria lead to a dramatic increase in 
rejected targets for |6| < 30, which is what motivates Criterion (0). 

Criterion (I) restricts the sample to stars with modest w-band photometric errors. Cri- 
terion (II) eliminates many stars, while preserving most HVS. Given the 4masyr _1 errors, 
criterion (III) eliminates a large fraction of stars surviving (II) while only eliminating 7% of 
HVS. Criterion (IV) ensures a minimum velocity in a direction consistent with an origin at 
the GC. Criterion (V) is the most important. It eliminates the overwhelming majority of 
halo stars and selects for metal-rich stars. The color index (u — g)o — 2.5(g — r)o "typically" 
(i.e., at g = 19.4) has errors of 0.105 mag, which means that if the color-color relation (see 
Fig [T|) is within a few hundredths of a magnitude of the threshold, then roughly half of the 
stars will survive this cut. Because of the relatively large error in this index, it is essential to 
keep the boundary away from the bulk of the metal-poor population, if one wants to avoid 
massive contamination. Hence, the combination of large errors and requirement of avoiding 
contamination dictates a color-color cut that will eliminate about half the HVS. 

We individually inspected a 20% random sample of the HVS targets rejected by the 
SEGUE-2 proper-motion selection algorithm on digitized POSS-I and POSS-II plates. We 
found that essentially all of those rejected because they had less than 4 astrometric epochs 
were, in fact, spurious targets: they had nearby (< 5") neighbors that were blended in 
POSS-I, thus corrupting the proper-motion measurements. Of the rejected targets that 
passed this test, almost none had any recognizable problem. These constitute 11% of all 
targets (excluding those rejected for too- few epochs). In addition, we estimate that < 2% 
of HVS (i.e., random field stars) would have a close neighbor that would lead to rejection 
as above. Hence, we apply a 0.87 (= 1 — 0.11 — 0.02) correction factor to our completeness, 
due to SEGUE-2 proper-motion selection. 

Our color selection is motivated by the requirement that we use minim al fibers per 



plate , so as not to interfere with any other program (of which there are many; lYanny et al. 



20091 ). To assess the quality of our p hotometric pre-election we rely on guidance from 



the theoretical isochrones computed by lAn et al.l (120091 ) . Figure [T] shows several tracks of 



stars plotted in color-magnitude and color-color space at several ages for a solar metallicity 
isochrone and a super-solar isochrone. As can be seen from the figure, our selection recovers 
solar metallicity objects very well for colors 0.4 < (g — r) < 0.5. At redder colors, the 
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isochrones formally fall below our selection, but by an amount that is smaller than the 
color-index error. At bluer colors, we retain substantial sensitivity only to stars above the 
main-sequence turnoff. For higher-metallicity stars, the isochrones shift upward in the color- 
index panels, improving our sensitivity at all colors. We are therefore sensitive to solar 
and super-solar metallicity ejectees. These curves shift downward at lower metallicity (not 
shown) away from our color-index selection region, and we lose sensitivity to sub-solar stars, 
except above the turnoff. This will be important to consider when computing our sensitivity 
to ejectees in Section [31 At solar metallicity, the isochrones straddle the color-index threshold 
in the region 0.4 < (g — r)o < 0.6, deviating by an amount that is small compared to the 
0.105 color-index threshold over the entire range. Hence, we obtain roughly 50% sensitivity 
to solar-metallicity G stars, at a cost of only about 1 fiber per plate, thus not substantially 
interfering with any other program. And, at higher metallicity, the recovered fraction is 
substantially higher. 



3. Results 

3.1. Where are the Old-Population HVS? 

The selection described above resulted in a total of 361 target stars within 181 SEGUE- 
2 fields. The velocity distribution of these targets is shown in Figure [2j As can be seen 
from the velocity distribution, while there are high-velocity objects, none of the targets is 
moving in excess of, nor even close to, Galactic escape speed (e.g., Xue et al. 2007). The 
selection therefore uncovered no metal-rich ejectees from this carefully constructed sample. 
Our null result may reflect a true dearth of old-population HVS in the halo, or it may reflect 
a target selection algorithm that is too stringent. In order to determine which is the case, 
we compute our sensitivity below. 



3.2. Sensitivity to Metal-Rich Ejectees 

We begin with a simple order-of-magnitude estimate of our survey sensitivity. We first 
define what is meant by "sensitivity" . Consider a Galactocentric spherical shell of thickness 
Ar. A star moving with speed v out of the GC will spend a time At = (1000kms _1 /f) x 
(Ar/lkpc) Myr traversing the shell. The total sensitivity to objects at a given velocity is 
therefore the integral of this quantity over the survey magnitude range and volume probed 
including selection effects. For this initial example, consider targeting every star of a single 
stellar type, with absolute magnitude M g = 4.5, to the survey limiting magnitude of g = 
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20 over the full area of the SEGUE-2 survey. Such stars can be probed to distances of 
approximately d = 10 a2 ( s-Mfl+5 )pc = 12.5 kpc from the Sun. Each SEGUE-2 spectroscopic 
plate has an area of VL p i ate = 7 deg 2 , and there are a total of N p i ates = 181. The survey volume, 
as seen from the Sun, is therefore V pro b e d,e = N p i ate £l plate d 3 / '3 = 265 kpc 3 . This corresponds 
to a Galactocentric shell thickness of Arcc = V P robed,e/ '^d 2 g1 where d g is the corresponding 
Galactocentric distance of the survey limit. At a fixed velocity, we can compute the time 
an ejectee spends in this shell, and therefore, the sensitivity as defined above. For a fiducial 
survey direction in which d g corresponds to 10 kpc, and a fiducial velocity of 500kms _1 , the 
star spends a total of 0.42 Myr traversing the shell, and the survey would be sensitive to 
stars of this type if they were ejected at a rate of ~ 2.3 Myr -1 . 

To compute the true survey sensitivity, we begin by fixing the absolute magnitude 
at a range of values (as above), and also consider a continuous range of ejection energies 
(parameterized by the HVS velocity at the time it passes the solar circle). For each of the 
SEGUE-2 plates and for each apparent magnitude in our selection interval, we compute d, d g , 
and the local velocity, and then simulate measurement of the star, allowing for proper-motion 
errors. The resulting "measurements" are then fed into the same algorithm that was used 
to select HVS candidates. Summation over all fields and apparent-magnitude intervals then 
yields the total sensitivity (in Myr). We no te the sharp contrast i n character between our 



completeness corrections and those made by lKollmeier et al.l ( 120091 ). That sample, although 



huge in an absolute sense, contained only a small fraction of potential HVS stars because it 
was not pre-selected. Hence there was a large completeness penalty. The present sample has 
stringent selection, so that the only completeness factors are those due to photometric errors 
(which scatter HVS in and out of our selection box) and the SEGUE-2 proper-motion-error 
completeness factor of 0.87. Hence, our completeness factor (computed rigorously within our 
code) is of order 50%. The result is shown in Figure [3j At M g = 4.5 and v Q = 500 km s -1 , 
the sensitivity is 0.16 Myr, consistent (within a factor of 3) with our order of magnitude 
calculation. 

Finally, to determine the sensitivity of the survey, we integrate the functions shown in 



Figure [3] over isochrones ( lAn et al.l 120091 ) of various ages and metallicities. That is, at each 
mass step (covering Alogm), with specified g — r color and ugr color index, we multiply 
A log m by the value indicated in Figure Eland by the fraction of stars that survive our color 
selection (given the photometric errors), and sum over the entire isochrone. This yields a 
sensitivity S (in Myr-dex). We therefore compute our sensitivity for isochrones of fixed age 
and a range of mass. Since we detected no HVS in our survey, we have 1 — e~ n n°/n\ — » 95% 
confidence that there are fewer than n = 3 expected detections over the range of mass and 
age to which we are sensitive. That is, we obtain a lower limit T > 3S 1-1 , where T is the rate 
of HVS ejections per Myr per dex. In Figure H] we plot this lower limit on F as a function of 
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population age for a range of metallicities at a fixed velocity at the solar circle of 500 km s 1 
As can be seen from the figure, our strongest limits apply to super-solar metallic ity stars with 



ages b etween 3-5 Gyr. For comparison, we also show the rate derived from the iBrown et al. 



(120091 ) survey, which is sensitive to stars of ages of roughly 150 Myr. Rather than limits, 
that survey has secure detections, so we plot it as a single point in the diagram. However, 
it is clear from this figure that our survey probes a very different set of stellar parameters 
relative to Brown et al. (2009). 



3.3. Limits on Metal-Rich Ejectees 

Our survey uncovered no metal-rich hypervelocity stars. As we have just demonstrated, 
our well-defined selection allows us to convert this null detection into an upper limit on the 
ejection rate as a function of stellar type. In Figure [5j we show our ejection limits as a 
function of velocity at the solar circle for several population ages and metallicities. We also 
show for comparison the results of the Brown et al. (2009) survey over the velocity range 
these HVS were observed. At sub-solar metallicities, our limits are no more constraining 
than the K09 results (slightly less stringent in fact; K09 found that there can be at most 
100 times the ejection of old-population stars as young population stars and we find a limit 
of 120). However, at solar metallicity and above, our limits are significantly stronger. For 
these stars, there can be no more than a factor of 30 times more solar-metallicity, low-mass 
stars relative to B-star ejectees. 

At ~ 500 km s -1 our limit on the ejection of solar-metallicity stars from Figure |5] is 
roughly r 5Gyr ' solar < 410(Myr - dex) -1 . At [Fe/H]= +0.2, our limits are improved to 

r 5Gyr,super-solar < 1 75 (Myr - dex)~\ 



3.4. Comparison to Young HVS 



The IBrown et al.l (120091 ) sample has 14 detected HVS objects over a total survey area 



Af2 = 5000 deg 2 , which corresponds to an ejection rate of roughly 1.75 Myr -1 over a mass 
range of roughly 3-4M (0.12 dex), or r 100Myr = 14 (Myr — dex) -1 . For a Salpeter mass 
function with shape dN/dlogM oc M -1,35 , we naively expect a factor of 5 more stars per 
dex at IMq compared to 3.5Af available at the GC for ejection. Of critical importance, 
however, is the relative life times of these two populations. Our sensitivity is maximized 
to stars of ages 5 Gyr. The IBrown et al.l (120091 ) survey is targeted for stars with lifetimes 
of roughly 150Myr. Therefore we must consider this additional factor of 33 times more 
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stars available for ejection, were they in fact accumulating at the GC for 5 Gyr. Therefore, 
provided the mass function is roughly Salpeter at the GC and that stars were accumulating 
for at least 5 Gyr, there should be ~165 times more old-population stars (which our survey 
targets) relative to B-stars available for ejection. Our limits place the ejection rate at no 
more than ~ 30 times more old-population ejectees relative to the known population of 
young ejectees. These stars are either not present at the GC, they have sub-solar metallicity, 
or the ejection mechanism quite strongly prefers high-mass stars at high velocity, by a factor 
of roughly 5.5. 



4. Discussion and Conclusions 



We have implemented a well-defined target selection algorithm within the SEGUE-2 
survey to search for metal-rich old-population hypervelocity stars. We targeted over 300 
stars and found no hypervelocity stars. This null result allows us to place new limits on the 
ejection of metal-rich old-population HVS from the GC as a function of mass, metallicity, 
velocity, and population age. Our results imply that either the mass function at the GC 
is top-heavy, or the mechanism for ejection of HVS significantly favors massive stars at the 
level of 1:5.5, or that the stars being ejected from the GC are sub-solar metallicity. 

It is certainly plausible that the ejection mechanism favors high-mass stars, or more 
precisely, ejects low-mass stars at velocities too low to be detected by our survey. It is 
possible that low-mass stars are being ejected at v elocities below our survey threshold of 



400 km s , and some mechanisms, e.g., iHillsl ( 19881 ). would predict t hat low mass b i naries 



would be ejected at lower velocity. Indeed, detailed calculations by iKenyon et al.l (120081 ) 
show that at 1M Q , the predicted velocity distribution at 400 km s -1 is strongly suppressed 
relative to more massive stars - our result in Figure [2] can be quantitatively compared 
to their predicted radial velocity distribution to determine whether the observations are 
consistent with the predictions for bo und ejectees. Furthermore, hig h-mass stars are thought 
to have a binary fraction near unity ( IPinsonneault & Stanekl 120061 ). The binary fraction is 
i ndeed observed to be substantially less than unity at lower mass in the solar neighborhood 
( jDuquennoy fc Mayor! Il99ll ; lLadal 120061 ). and it is possible that conditions at the GC result 
in a similar F/G binary fraction. Therefore, the classic binary disruption mechanism could 
very well favor high-mass stars in this scenario. Relaxing the assumption of a universal mass 
function, however, it i s also p ossible that there are simply very few old stars at the GC. For 
example, iPerets et al.l ( 120071 ) demonstrate that, for reasonable assumptions about the mass 
function and old-star b inary-fraction, the ejec tion rate of 1M Q stars should be approximately 



5 x 10 7 yr 1 (but see lYu &: Tremaind 120031 ) . and could be substantially enhanced by the 
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presence of a secondary massi ve per t urber as envisioned by e.g., iPolnarev fc Reed (11994') . 
Baumgardt et al.l (120061 ) . and iLevinl (120061 ) and not included in the iKenyon et al.l (120081 ) 
calculations. 

One mechanism that has been suggested is that the HVS have nothing at all to do with 
the black hole at the GC and the dynamics there, b ut rather they reflect debris of a satellite 
with pericenter close to the GC ( lAbadi et al.l 120091 ) and a large population of young stars. 
The anisotropy of the known HVS potentially supports this interpretation. A satellite galaxy 
would have the normal complement of old-population stars compared to young stars, and 
that would show up in our limits. While we do not see these stars, if the satellite system had 
low metallicity similar to the current satellites of the Milky Way, it would not be detectable 
in our survey. It is also possible that our SEGUE-2 pointings were not sufficiently aligned 
with the proposed direction of this debris. 

The new limits derived here, in conjunction with the known HVS detections at higher 
mass, provide an important constraint on models of HVS ejection. The target strategy 
implemented here represents an efficient and economical way to probe low mass populations 
ejected from the GC at modest spectroscopic follow-up cost. Future spectroscopic surveys 
with similarly stringent target selection could strengthen these limits. Deciphering the true 
distribution of GC ejectees as a function of mass, age, velocity, and position will continue to 
provide valuable insight into the still obscure Galactic center. 
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Fig. 1. — Color-magnitude and color-index diagram for theoretical isochrones computed 
by An et al. (2009) as a function of age and metallicity. Left panels show results for a 
[Fe/H]=0.0 isochrone; right panels show a [Fe/H]=+0.2 isochrone. Colors green, blue, and 
cyan show ages of isochrones corresponding to 1.2, 3.5, and 5.0 Gyr, respectively. Our color 
selection boxes are superposed on the isochrones in the lower panels as black rectangles. Our 
survey has sensitivity to solar and super-solar metallicity stars. 
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Fig. 2. — Galactocentric velocity distribution of HVS targets from SEGUE-2. 
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Fig. 3. — Sensitivity to HVS as a function of velocity at the solar circle and absolute mag- 
nitude of stellar population. Curves from top to bottom show our sensitivity as a function 
of absolute magnitude ranging from M g = 4.4 (top) to M g = 5.0 (bottom) in increments of 
0.1 magnitudes. 
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Fig. 4. — Limits on the ejection of old-population HVS as a function of metallicity and age, 
for velocities at the solar circle v & = 500 km s -1 . Solid, dotted, short-dashed, long-dashed, 
and dot-dashed lines correspond to metallicities of [Fe/H]= +0.2, +0.1, 0.0, —0.1, and —0.2, 
respectively. For comparison, we show the Brown et al. (2009) survey results, which probe 
a very different age regime compared to our survey. Our ejection limits are most stringent 
for metal-rich stars; our sensitivity to sub-solar metallicity is significantly degraded. 
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Fig. 5. — Limits on the ejection of old-population HVS as a function of metallicity, age, 
and velocity. Solid, dotted, dashed lines correspond to metallicities of [Fe/H]=0.0, —0.2 and 
+0.2, respectively. Colored lines show different population ages of green, blue, cyan to 1.2, 
2.5 and 5.0 Gyr, respectively. Each point of the solid curves is formed by integrating over a 
set of corresponding points at the same velocity in Figure [31 weighting by the fraction that 
survive the color selection for each isochrone, and dividing 3 by the result. For comparison, 
we show the Brown et al. (2009) survey results. 



